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Atomic quantum superposition state generation via optical probing 
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We analyze the performance of a protocol to prepare an atomic ensemble in a superposition of two 
macroscopically distinguishable states. The protocol relies on conditional measurements performed 
on a light field, which interacts with the atoms inside an optical cavity prior to detection, and we 
investigate cavity enhanced probing with continuous beams of both coherent and squeezed light. 
The stochastic master equations used in the analysis are expressed in terms of the Hamiltonian of 
the probed system and the interaction between the probed system and the probe field and are thus 
quite generally applicable. 

PACS numbers: 42.50.Dv, 42.50.Pq 



I. INTRODUCTION 



The ability to prepare quantum mechanical systems 
in superposition states is important both in fundamental 
studies of quantum mechanics and in technological appli- 
cations. Generation of such superposition states of light 
by conditioning on the outcome of a measurement has 
been achieved experimentally [l|, H, but light fields 
are difficult to store for long periods of time, and it is 
desirable to be able to prepare trapped atoms in super- 
position states as well. First steps in this direction have 
demonstrated generation of superposition states of up to 
6 trapped ions [1] and W states of up to 8 trapped ions [f| 
through interactions with suitably chosen light pulses. In 
the present paper we consider trapped atomic ensembles 
and investigate the possibility to prepare superposition 
states by means of optical quantum non-demolition mea- 
surements. 

The idea to probe the state of an atomic system by 
allowing the system to interact with a light field and 
then subject the light field to measurements is very use- 
ful. It has, for instance, been used to observe quantum 
jumps of single ions [H, Another application is the 
generation of spin squeezed states by means of quantum 
non-demolition measurements @, and related proposals 
use spin squeezing to improve the precision of atomic 
clocks [13, fll[ and mag netometers lip . Several recent 
experiments [l3l 0, [TJ EE G3, 0, E3 have focused on 
the possibility to trap cold atoms and Bose-Einstein con- 
densates inside high-finesse optical cavities, because the 
cavity enhances the light-atom interaction strength, and 
the strong-coupling regime, where the coherent dynamics 
takes place on a faster time scale than the dissipative dy- 
namics, has been reached. In addition, optical measure- 
ments can be used to investigate the atom statistics of 
atom laser beams (2fjj and the atomic population statis- 
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tics of optical lattices [2l|, |2J, \2M ■ 

In the following, we demonstrate that for the strong- 
coupling parameters obtained in [l7| , a measurement pro- 
cedure similar to the quantum non-demolition measure- 
ment used to generate spin squeezed states can be used 
to g enerate quantum superposition states, as outlined in 
|24| (see also [25[ for a related proposal). In Sec. HT1 we 
explain the state preparation protocol, and we provide 
the time evolution of the state of the atoms and the light 
field. The performance of the protocol is analyzed in 

Sec. nu 

It is convenient to use a continuum coherent state from 
a laser as the probe field, but it is also interesting to inves- 
tigate the additional possibilities that arise, if we choose 
to probe the system with different quantum mechani- 
cal states of light. An example of free-space probing of 
an atomic ensemble with single-mode photon subtracted 
squeezed vacuum states has been investigated in [2fj| , and 
in Sec. IIVI of the present paper we provide general tools 
to analyze probing with a continuous beam of squeezed 
light. Despite the infinite dimensionality of the Hilbcrt 
space of the light field it is possible, for the case at hand, 
to rewrite the equation governing the time evolution of 
the state of the atoms and the light field into a finite set 
of coupled stochastic differential equations, and this al- 
lows us to obtain numerical results for the full dynamics. 
Section IVl concludes the paper. 



II. PROBING PROCEDURE 

To probe atoms inside an optical cavity, we send a 
probe beam, which is initially in a continuum coherent 
state, into the cavity and observe the light reflected from 
the cavity with a homodyne detector as shown in Fig. [T] 
We consider atoms with two ground state levels, which 
we treat as a spin- 1/2 system, and we assume that the 
cavity field couples one of these levels non-resonantly to 
an excited state. Assuming that all atoms couple identi- 
cally to the cavity field, which may be achieved e.g. with 
a Bose-Einstein condensate, and eliminating the excited 
state adiabatically, the Hamiltonian for the light-atom 
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FIG. 1: The atoms are placed inside an optical cavity and are 
probed with a light beam, which is in a continuum coherent 
state before it enters into the cavity. The probe light reflected 
from the cavity is observed with a homodyne detector and the 
atomic quantum superposition state is generated for certain 
possible measurement readouts. 
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FIG. 2: Real and imaginary parts of a n for g — 2n x 215 MHz, 
A = 2tt x 10 GHz, ki = k = 2tt x 106 MHz, 4ki/3 2 /« 2 = 0.01, 
and n = —100, —99, . . ., 100. The points lie on a circle with 
center (^/kI/3/«, 0) and radius ^/kIP/k. 



interaction takes the form 



H = hgaJaJ. 



(1) 



where g = g 2 /A, g is the single-atom coupling strength 
on the optical transition, A = lu — ujq is the detuning, uj 
is the frequency of the light field, loq is the frequency of 
the atomic transition, a is the field annihilation operator 
of the cavity field, J z is the z-component of the collective 
atomic spin vector J = J^ji, and ji is the spin vector 
of the zth atom. We have here ignored the possibility 
of spontaneous emission from the excited state, which is 
valid in the strong coupling limit. We return to this point 
in Sec. UTT. 

To understand why superpositions of very different 
quantum states are obtained for certain measurement 
outcomes, if the phase of the local oscillator is chosen 
appropriately, we first consider the case, where the ini- 
tial state of the atoms is a single eigenstate of J z with 
eigenvalue hn. In this case the unitary time evolution op- 
erator corresponding to the above Hamiltonian reduces 
to a phase shift operator, and provided the cavity field 
is in the vacuum state at time t = 0, one finds that the 
cavity field at time t is a coherent state with amplitude 



a n (t) 



>/2+*»#)(*-*') /3( t ')dt' ', (2) 



where K\ is the cavity decay rate due to the input mirror, 
k = K\ + Kioss.i is the total cavity decay rate including 
loss, and (3(t') is the amplitude of the probe beam, i.e., 
\P(t)\ 2 At is the expectation value of the number of pho- 
tons in the beam segment of length cAt, which arrives 
at the input mirror at time t. Note that ([2]) is indepen- 
dent of the actual outcome of the homodyne detection, 
which appears because the quadrature noise of the co- 
herent state leaving the cavity is uncorrelated with the 



quadrature noise of the coherent state inside the cav- 
ity. Choosing Pit') to be real, it follows from $2$i that 
a n (t) — a*_ n {t), and a measurement of the x-quadrature 
of the cavity field is thus unable to distinguish between 
n and — n. This suggest that for a general initial state of 
the atoms, the homodyne detection slowly projects the 
atomic state onto a superposition of states with the same 
value of |n|. 

If we choose (3{f) to be time independent (in a frame 
rotating with the relevant bare cavity resonance fre- 
quency), a n (t) assumes the steady state value 



2y/K\j3 1 — ling I k 
k 1 + An 2 g 2 1 k? 



(3) 



after a transient time of order a few The ampli- 

tude a n is plotted in Fig. [2] for the values of g and k 
in Ref. [17[ and different values of n. For small 2\n\g/K, 
a n is simply rotated by an angle proportional to n com- 
pared to o;o, but when the rotation angle increases, the 
amplitude of the cavity field drops, because the probe 
field is no longer on resonance due to the phase shift in- 
troduced by the atoms, and as 2\n\g/K tends to infinity, 
a n approaches the origin. Thus, \a n — a_ n | is small for 
2\n\gJ K 3> 1, which means that the rate of decoherence 
between |n) and | — n) due to cavity losses and detec- 
tion inefficiency is reduced. To illustrate this explicitly 
we note that if the initial state of the atoms is an equal 
and pure superposition of n and —n, then it follows from 
Eq. iJSJ) below that the purity of the atomic state evolves 
as 



1 



1 + exp 



4mP 2 n 2 g 2 



[{k/2) 2 + n 2 g 2 } 2 

X[(K-TIKl)t+l]\\ (4) 
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where Pa±(t) is the atomic density operator obtained by 
tracing out the cavity field, r\ is the detection efficiency, 
and we have neglected transient terms proportional to 
exp(— Kt/2) and to exp(— nt) in the exponent. In ad- 
dition, the smaller number of photons in the cavity for 
large 2\n\g/n leads to a decrease in the rate of sponta- 
neous emission from the excited atomic state. 

The stochastic master equation for the setup in Fig. [1] 
was derived in Ref. [H| for ngr < 1, w < 1, and 
Y^kTI/JIt <C 1, where r is the round trip time of light 
in the cavity, and these requirements are all fulfilled for 
the parameters in Fig. [5] and a cavity length of 38.6 pin, 



as in Ref. [17]. In the following, we consider the case, 
where the initial state of the atoms is a coherent spin 
state pointing in the x-direction, and for notational sim- 
plicity, we thus restrict ourselves to a basis consisting 
of simultaneous eigenstates of J 2 and J z with total spin 
quantum number J = N/2, where N is the number of 
atoms. This is possible because the Hamiltonian com- 
mutes with J 2 . If we also assume j3 = /3* and a local 
oscillator phase, which corresponds to a measurement of 
the x-quadrature of the cavity field, the solution takes 
the form 



J J 



p(t)=Y, Y\ Cnm(0)exp -£/ (\a n {t')\ 2 + \a m (t')\ 2 -2a n (t')a* m (t'))dt' + ^ f (o»(t') + a* m (t')) dy' 

n=~ Jm =-J I Z JO JO 



- I (a n (t') + a* m (t')f dt' - i^T / 0{t') (a n (t r ) - a* n (t') - a m {t') + a* m (t')) dt' 



t (t)){a m {t)\®\n){m\/ J2 c ii(°) ex P 

1 q=-J 



VKl / {a q (t') + a*(t'))dy' - -rjm / (a g {t') + a* q {t')fdt 



> (5) 



where p(t) is the density operator describing the state 
of the cavity field and the atoms, \a n (t)} is a coherent 
state with amplitude @, \n) is the eigenstate of J z with 
eigenvalue tin, and dy' is proportional to the photocur- 
rent measured in the time interval from t' to t 1 + dt' (see 
Ref. for details). The coefficients C m „(0), normal- 
ized according to ^] C w (0) = 1, are determined from 
the initial state of the atoms, and for a coherent spin 
state pointing in the x-direction 



C ram (0) 



(2J)! 



4 J y/(J + n)\(J - n)\{J + m)\{J - m)\' 



(6) 



We note that the trace over the atomic state of p(i) is an 
incoherent sum of the states |a ra )(a ra |, and if the differ- 
ent coherent states are given appropriate weights, Fig. [2] 
may be regarded as a phase space representation of the 
cavity field. If, instead, we trace out the light field to 
obtain p a t{t), we observe that the weights of |«)(n| and 
| — n)(— n\ are equal at all times if they are equal at the 
initial time, as expected. Note also that the time evolu- 
tion of the diagonal elements of p a t{t) is independent of 
the off-diagonal elements of p a t (t) . 

The density operator at time t depends on the results 
of the continuous measurement, and the probability den- 
sity to observe a given sequence dy±, dj/2, dyu of 



measurement outcomes in the interval from to t is 



P{dy 1 ,dy 2 , . . -,dy M ) = Y[ 



k=l 



\ 1 ( 


dyl\] 


, exp 

_V2irdt V 


2dt ) 



M 



x C OT (0)exp }2^qi^Y, Bje l a i(( k - 1)#)]<*M 

q=-J ^ k=l 



t 

-2<r)K X I [Re(a q (t'))] 2 dt' 
lo 



(7) 



where dt = t/M and M is very large. One could now 
use Eq. ([7]) to derive the joint probability density for 
the 2{N + 1) stochastic variables J Re(a n (t'))dy' and 
f Q Im(a n (t'))dy' , but for constant (3 and a total probing 
time, which is large compared to k , we may approx- 
imate a n (t') by its steady state value a n . In this case, 
p(t) only depends on the stochastic variable Y = f Q dy', 
which has probability density 



pp) = E 

q=-J 



c qq (Q) 

f2^i 



exp 



(Y - 2^qK{Re{a q )tf 



2t 



(8) 



III. PERFORMANCE OF THE STATE 
PREPARATION PROTOCOL 

The ability to distinguish the atomic states \n) and 
\m) is determined by the distance |Re(a„ — a m )\ be- 



4 



tween the corresponding coherent states measured along 
the x-axis in phase space. A glance at Fig. [2] thus re- 
veals that we may conditionally produce states, for which 
the distribution (n|p a t(t)|n) consists of two peaks sepa- 
rated by a region with (n|p a t(t)|n) ~ 0, after a relatively 
short probing time, while it takes significantly longer to 
produce a superposition of states with a single value of 
\n\. To estimate the required probing time to distin- 
guish between |0) and | ± J), we note that the probabil- 
ity density ([5]) is a sum of 2 J + 1 Gaussians with mean 

For 



values 2^/i]KiKe(a q )t and standard deviations \/t. 
Ag 2 J 2 /K 2 > 1, Re(a±j) » (see Eq. (J3])), and the Gaus- 
sians corresponding to n = and n = ± J thus begin to 



separate after a probing time 



to 



r]Ki 4ki/3 2 



(9) 



which evaluates to 150 ns for the parameters in Fig.[^]and 
i] = 1. To clearly distinguish between n = and n = ±J, 
we need to choose a probing time, which is somewhat 
larger than this, and Fig. [3] suggests that t « 1CP 6 s is 
a reasonable choice for the considered parameters. Note 
that t qs is independent of the number of atoms as long as 
4g 2 J 2 /n 2 3> 1. Note also that t sa 10~ 6 s is large com- 



pared to k = 1.5 x 10 1 
state assumption Kt ^> 1. 



s, which justifies the steady 
The maximal probing time 



J 



is limited by decay processes, and, in particular, we re- 
quire that t is small compared to the average time t sp 
between spontaneous emission events, which may be ap- 
proximated as the inverse of the product of the atomic 
decay rate T, the probability for a single atom to be in 
the excited state, and the number of atoms in the state, 
which interacts with the light field, i.e., 



l + 



N 



-i -i 



A 2 + (r/2) 2 2 

\ (10) 

where h is a typical value of n. The smallest approximate 

= and is 6 x 10~ 5 s for the 
2ttx6MHz (asinRef. [13]), 



value of t sp is obtained for h 
parameters in Fig. [2] and T = 



while wc obtain t sp ps 2 x 10 s if we choose |n| to be 
the most probable value of \n\ after a probing time of 
t = 10~ 6 s and a measurement outcome Y — s 1 / 2 (see 
Eq. (| 1 5(1 below). It is thus realistic to finish the probing 
before one of the atoms decay by spontaneous emission, 
and we hence neglect atomic decay throughout. Note 
that this conclusion is independent of j3 since the rate of 
gain of information and the rate of spontaneous emission 
events are both proportional to the probe beam intensity. 

An important measure of the quality of the generated 
atomic states is the purity 



Tr (/&(*))= £ £ |C nm (0)| 2 exp 

ri — — J m— — J 

+ Vwl I (a n (t') + a* n (t') + a m (t') + a* m (t'))dy' 
Jo 



Mt) " "mWI 2 - (« - rjKi) I \a n (t') - a m (t')\ 2 dt' 



rim 



(K(t') + a* n (t')f + Mt') + a* m (t')) 2 )dt 



q=-J 



/rj«r / (a g (t') + a* q (t'))dy' 
Jo 



(a,(f) + a q {t')) 2 dt 



(11) 



r 



The factor exp(— \a n (t) — a m (t)\ 2 ) appears, because the 
cavity field is traced out to obtain the atomic state, but 
this reduction in purity may be eliminated by turning off 
the probe beam after the desired probing time, whereby 
the cavity field quickly decays to the vacuum state and 
the density operator factorizes. If, in addition, the ini- 
tial state is pure, and there are no losses (t]Ki = k), the 
atoms are seen to end up in a pure state, as expected. 
The purity is plotted as a function of the integrated photo 
current Y for r\K\ — 0.9 K in Fig. 01 neglecting the tran- 
sient of a n (t) , and the relevant range of Y- values may be 
read off from Fig. [5] For large positive values of Y, the 
purity is close to unity, because the atomic state |0)(0| 
has a large weight. The purity is smaller for intermedi- 
ate values of Y, because states with Re(a„) w Re(«o)/2 
have large weights, and the distance in phase space be- 



tween a n and a_„ is large for these states. For low values 
of Y the purity increases, because the states with large 
weights are closer to the origin, and |a n — a_ n | is thus 
smaller. The purity for low Y is different for different 
numbers of atoms, because \aj — a_j| depends on J. If 
the probe field is turned off at some point to ensure that 
a n (t) — at the final time, it follows from that the 
atomic state only depends on (3 through j3 2 dt and (idy 
(note that a n is proportional to /?), and a scaling of /3 
thus corresponds to a scaling of the probing time and the 
measured photo current. 

The distance d between the two peaks in Fig. [3jb) 
may be used as a measure of the difference between the 
two states in the generated quantum superposition and 
is thus a second parameter to characterize the quality of 
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FIG. 3: (Color online) Spin Q-function [23, ^ on the Bloch 
sphere (a) and diagonal of the atomic density operator (b) for 
the state produced at t = 10 -6 s, assuming Y = 5 x 1(F 4 s 1 ^ 2 , 
J = 50, rj — 1, and otherwise the same parameters as in Fig. [2] 
The crosses in (b) indicate the diagonal of the atomic density 
operator at the initial time. 



the states. For large N we may approximate ((6]) by 



C„ m (0) = \/ — exp 



n 2 + to 2 
N 



(12) 



and it thus follows from Eq. ([5]) that the largest value of 
{n\p a t(t)\n) is obtained for n — n p , where 



1 + 



N 



y/a (3V3abt + bVbY 3 + 27at 2 ) V3 - bY 



V3 (3V3abt + bVbY 3 + 27at 2 ) 

and we have defined 

, 2^/3 Ag 2 N 

a = ^JrjKi x — and o = ^ — . 



1/3 



(13) 
(14) 



FIG. 4: Purity of the atomic state as a function of the in- 
tegrated photo current Y for t — 10 -6 s, rjKi = 0.9 k, 
g = 2ir x 215 MHz, A = 2tt x 10 GHz, « = 2tt x 106 MHz, 
and iK!f3 2 /K 2 = 0.01. 
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FIG. 5: Probability density to measure a given value of Y 
for the same parameters as in Fig. U When the number of 
atoms increases the probability to observe lower values of Y 
increases because extra coherent states with small amplitude 
are added to the left in Fig.[2]and the width of the distribution 
(n|p a t(0)|n) at the initial time is broader. 



We plot n p in Fig. and it is seen that n p increases with 
decreasing Y, because small values of Y favor atomic 
states with |n| close to J. For Y — we obtain the 
particularly simple result 



Ag 2 n 2 p 



1 



2r/K 



A Kl f3 2 Ag 2 N 



1/3 



(15) 



and, in this case, d scales as TV 1 / 6 for large N, which 
means that the relative distance d/N between the peaks 
decreases with N. Note, however, that the probability to 
measure Y = is different for different values of J. The 
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FIG. 6: Value of \n\, which maximizes {n\p at (t)\n) , as a func- 
tion of Y for the same parameters as in Fig. [4] The points 
are obtained directly from Jjj), while the lines represent the 
approximation ([13}. 



ity field obtained for atomic states with n ss J and the 
cavity field obtained for n ss —J, in order to decrease 
the decoherence rate due to light field losses, and, at 
the same time, increase the difference between the cavity 
field obtained for atomic states with n rs and the cav- 
ity field obtained for n « ±J, in order to decrease the 
required probing time. The former may be accomplished 
by setting /3 = 0, but this choice also prevents us from 
retrieving any information about the state of the atoms 
from the measurement records. If, on the other hand, 
we squeeze the probe beam before it enters into the cav- 
ity containing the atoms, the cavity field is different for 
different values of n even if (3 — 0. This suggests that 
it might be possible to increase the purity of the final 
atomic states by using a squeezed vacuum probe field, 
and we thus investigate this possibility in the following. 
For the sake of generality, we shall, however, not put 
= until Sec. HVDl 



peak-to-peak distance d may be increased, by decreasing 
4^ 2 /k 2 , but since this moves the coherent states in Fig. [2] 
to the right around the edge of the circle, this will also 
decrease the final purity if losses are present. 

IV. PROBING WITH SQUEEZED LIGHT 

Considering the conclusions from the previous sections, 
we would like to increase the similarity between the cav- 



A. Equation of motion 



A setup to probe the atoms with squeezed light is de- 
picted in Fig. [7] and basically consists in passing the probe 
beam through a cavity containing a pumped nonlinear 
medium. Combining the theory developed in Refs. [24| 
and (29[, we derive the stochastic master equation 



p(t + dt) = p(t) - -[H, p{t)]dt - - [(e(c t ) 2 + e*c 2 ), p(t)] dt + - («i + Ki oss ,i) (-a)ap{t) - p{t)a) a + 2ap{t)tf) dt 

+ K ioss,2) (— $cp(t) — p{t)c) c + 2cp(t)c t ^J dt + iy / KiK2(ci , cp(t) — p(i)c^a + ap(t)c^ — cp(t)a))dt 

+ i^[~(3[a\ p(t)]dt + [a, p{t)]dt - ^(3[c* , p(t)]dt + ^7 2 f3* [c, p{t)]dt 

+ { - e-^V^T M*) - Tr P(M - e^V^T [ P {t)a) - Tr (p(t)fit) p ( t )] 

+ ie-^yfqi^ [cp(t) - Tr (cp{t)) p(t)} - ie^y/^ [p{t)$ - Tr {p{t)c)) p{t)] }dW. (16) 

I 



Here, pit) is the density operator representing the state 
of the atoms, the light field in cavity 1 (annihilation op- 
erator a), and the light field in cavity 2 (annihilation 
operator c). (We assume that each mirror gives rise to 
a phase shift of 7r/2 on the reflected field, and the phase 
of the light field thus assumes four different values along 
the length of each cavity. The operators a and c re- 
fer to the field passing through the atomic cloud and 
the crystal, respectively.) H is the Hamiltonian for the 



light-atom interaction and the internal dynamics of the 
atoms, e is the nonlinear gain coefficient in the crystal 
(the time evolution operator corresponding to one round 
trip in the cavity is U sq — exp{— i[e*c 2 + e(c' f ) 2 ]r2/2}, 
where ti is the round trip time of light in cavity 2, and 
l e l/( K 2 + Kioss,2) = 1/2 at threshold), k\ and k 2 are cavity 
decay rates due to the cavity input mirrors as indicated 
in the figure, ki OS s.i and ki OS s.2 are additional cavity de- 
cay rates due to cavity losses, j3 is the amplitude of the 
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FIG. 7: Probing with a squeezed state may be achieved by 
including an optical parametric oscillator in the setup. The 
input beam may be either a continuum coherent state or vac- 
uum, and the light field is squeezed when it passes through 
the pumped nonlinear crystal. 



coherent states along x that we use as initial states are 
not time-reversal eigenstates, but are eigenstates of the 
product of time-reversal and a rotation by ir around J z , 
exp(— iirJ z )T corresponding to a reflection in the xy- 
plane of the spin. Since T and exp(— i-n J Z )T are both 
symmetries of the propagation, our initial state evolves 
into a state with the same reflection symmetry between 
n and —n. 



C. Solution of the equation of motion 

For the Hamiltonian |T]) , which does not by itself give 
rise to transitions between the different |n) states, we 
may write the solution of (fT6|) as 



m = E E 



(19) 



n— — J rn=—J 



where p nm (t) are operators acting on light fields in both 
cavities, i.e., on a two-mode system. In the following, we 
use a mixed Wigner function density operator represen- 
tation and define W nm (xx,p±, X2,P2, t) as the (two-mode) 
Wigner transform of p n m(t): 



input beam, 4> is the phase of the local oscillator, and 
dW is a stochastic variable with a Gaussian probability 
distribution with mean value zero and variance dt. 



B. Symmetries 

In order to produce symmetric superposition states of 
the atoms, we should be able to choose the parameters 
such that our detection gives information on |n|, but not 
on the sign of n. Although intuitively clear from the 
physical setup, to formally prove that this is possible, 
we consider a variation of the time-reversal operator [30| : 
Let T be an anti-unitary operator which acts on the atom 
and light field operators as follows [34j : 



W nm (xi,pi,X2,P2, t) = 



fjf' 1 = —j, fdf- 1 = a, fcf- 1 = 



(17) 



Writing dW = £dt and p(t+dt) = p(t) + £z[p(t)]dt, where 
the subscript £ emphasizes the explicit dependence on the 
random measurement outcome, we then have 

fp(t + dt)^ 1 = fp^f- 1 + C$ \f pi^f- 1 ] dt (18) 

provided that <fi — kn, k 6 TL (corresponding to detec- 
tion of the X\ quadrature) and that (3 and e are both 
purely imaginary (note the phase shift imposed on the 
input beam due to the presence of the squeezing cavity) . 
Because of (|18j) , time- reversal symmetry is preserved dur- 
ing the propagation for each realization of the noise. In 
addition, H of Eq. (TTJ) is obviously invariant under arbi- 
trary rotations around the J z spin axis. Now, the spin 



1 

4^ 



Tr fe&r+irii)^ -(J7r— t'7i)a+(Cr+«Ci)e t ~^ r ~ i ^ & p nm (t)^j 
e -V2ir HXl + V2 i > nrPl -V2i( t X2 + V2 i ( r p 2dVrdmd £ rd ^_ ^ 

At the same time, we also need to translate the light field 
operators in Eq. (116[) into operators acting on Wigner 
transforms (see [31| chapter 4). For a coherent input 
beam 



W nm (y,t) 



N nm (t) 



ir 2 ^/det(V nm (t)) 
cxp [ - (y - y nm (t)) T {V nm (t))- l (y - y nm (t))] , (21) 

where y = (xi,pi,X2,P2) T is a column vector of the 
quadrature variables of the field in cavity 1 and cav- 
ity 2, and, for each n and to, V nm (t) — V nm (t) T — 
V^ n {t) are symmetric and complex four-by-four matri- 
ces, y nm {t) = Vmnit) are complex four-by-one vectors, 
and N nm (t) = N^ n (t) are complex scalars. The problem 
thus reduces to a finite set of coupled stochastic differ- 
ential equations for V nm (t), y n m{t), and N nm (t), which 
may be solved numerically for given realizations of the 
noise dW . In addition, for imaginary e and /3 and real 
e l< ^ the time-reversal symmetry described above implies 
that the four-by-four matrices V— m< — n (t) and V nm (t) are 
related by a change of sign of the elements in the second 
and third rows followed by a change of sign of the ele- 
ments in the second and third columns, that j/_ m) _„(t) 
is obtained from y nm (t) by changing the sign of the sec- 
ond and third rows, and that N- m ^- n {t) — N nm {t) if 
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iV_ mi _„(0) = N nm (0). As expected, it thus follows that 
the Wigner function phase space representation of the 
field in cavity 1 is symmetric under reflection in the aq- 
axis and that the atomic states |n)(n| and | — n) (— n\ have 
equal weights at all times if they have equal weights at 
the initial time. 

The differential equations reveal that the evolution of 
V nm (t) is deterministic and that dV nm (t) / dt only de- 
pends on n, m, V nm {t), and the physical parameters 
of the setup, except (3. Numerically we observe that 
Vnm(t) quickly approaches the steady state solution, 
and it is thus possible to reduce the number of differ- 
ential equations significantly if the transient is negligi- 
ble. In contrast to the case without squeezing, y nm (t) 
depends on the measurement outcome for e ^ be- 
cause the noise of the field leaking out of the cavities 
is correlated with the noise of the field inside the cavi- 
ties, and the observations lead to a back-action on the 
mean value of the field amplitudes for each n as well 
as an update of the probabilities. In the unobserved 
case (77 = 0), Dnn{t) does, however, approach a steady 
state value, which for (j> = 7T, = i\(3\, and e = ilm(e) 
is ), \/2Im(a £ ,„), 0, -2|/3| V2Ka/(«2 + ^loss,2 t 

2Im(e))] T , where a t:n is given by the right hand side of 
Eq. ([3]), except that /3 is replaced by \ f3\ [2k 2 / (k2+ki oss 2 + 
2Im(e))-l]. 

The stochastic nature of y nn (t) prevents us from draw- 
ing a graph corresponding to Fig. [2] which is valid for 
all possible measurement outcomes, and instead we note 
that V nn (t) is the four-by-four covariance matrix of the 
two light field modes conditioned on the atomic state 
|n)(n| and plot the corresponding steady state error el- 
lipses of the cavity 1 field in Fig. [51 As in the unsqueezed 
case, extreme n values, \n\ ~ J, leads to cavity 1 being 
almost in the vacuum state because the effective reso- 
nance frequency of the cavity is shifted by the atoms and 
the probe beam is hardly coupled into the cavity at all. 
It is thus not surprising that the squeezing is very weak 
for these extreme n-values, as observed in the figure. 



D. Numerical results 

In Fig. [9] we plot numerical results for (3 = and var- 
ious values of the nonlinear coefficient e 35], and, for 
reference, we also plot results of probing with a coherent 
state field with 4ki|/3| 2 /k = 0.01. For nonzero e, the ob- 
tained state depends on the whole detection record and 
every realization of the noise dW leads to a unique state. 
It is thus more difficult to summarize the results. We 
have chosen to present a scatter plot of purity versus d/2, 
where ±d/2 are the positions in the peaks in the n dis- 
tribution of the atoms, cf. Fig.[3]Jb). Larger values of d/2 
thus correspond to a more "macroscopic" superposition 
state. 

For the strongest squeezing considered, e = ±0.05 !«2, 
we see that for all but the highest values of d/2, the pu- 
rity is lower than what is obtained with coherent state 




xx 

FIG. 8: (Color online) Steady state uncertainty ellipses of 
the light field in cavity 1 corresponding to the atomic state 
\n)(n\ for different values of n. The parameters are g — 2n x 
215 MHz, A = 2tt x 10 GHz, <f> = ir, e = ±0.2 IK2, Ki = n 2 = 
2tv x 106 MHz, Kioss.i = «:i oss ,2 = 0, and 77 = 0.9. The blue 
(dark) lines are for e = +0.2 iki and it is apparent that the 
state corresponding to n — is squeezed in the xi-quadrature 
and anti-squeezed in the pi-quadrature. As |n| increases, the 
squeezing is rapidly reduced and the uncertainty ellipse turns 
so that the minimal (maximal) uncertainty is no longer along 
xi (pi). The quadratures thus become correlated and back- 
action due to an x\ measurement also affects the pi mean 
value. The red (bright) lines are for e = —0.2 in2 and in that 
case, the n — state is anti-squeezed in the ^-quadrature and 
squeezed in the pi-quadrature. Note that the error ellipses are 
independent of and that, for the purpose of illustration, we 
have chosen a rather high value of |e|. 



probing like in Sec. IIIII We attribute this to the correla- 
tions between the field leaking out of the cavities and the 
field inside the cavities: In contrast to a coherent state, 
the purity of a squeezed state decreases if it is subjected 
to loss, and the stochastic movement of the mean values 
of the cavity field for different values of n may increase 
the distance in phase space between the cavity fields cor- 
responding to the same value of \n\. 

In fact, it is not quite fair to compare the e = ±0.05 IK2 
results to the coherent state results as for the chosen ob- 
servation time, t = 10 -6 s, the former has more peaked n 
distributions than the latter. This means that one could 
have reduced t in the squeezed case. Instead we keep t 
fixed and look at results for e = ±0.025 iK2, which for 
t = 10~ 6 s leads to n-distributions very similar to the 
coherent state results. As can be seen from Fig. [9j we 
can then preserve a high purity while still being sensi- 
tive enough to make the n distribution double-peaked. 
The highest purity is found for the phase of e where the 
n = light-state is squeezed in the p\ quadrature. A 
closer analysis reveals that for higher |n| the measure- 
ment back-action is then less likely to drive the light field 
amplitudes for ±n apart and thus the field leaking to the 
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FIG. 9: (Color online) Scatter plot of purity vs. d/2, where 
±d/2 are the positions of the peaks of the n-distribution of 
the atoms. The larger d/2, the more separated are the two 
components of the macroscopic superposition, cf. Fig. EJb). 
The results are for J — 50 and coupling and loss parame- 
ters like in Fig. [8] The input light and squeezing is on for 
t = 10~ 6 s, but we continue observations for an additional 
time of 15/ki ~ 2 x 10~ 8 s in order for the cavity light 
field to decay to vacuum. The symbols (black) indicate 
results for a coherent state probe with 4ki|/3| 2 /k 2 = 0.01. 
The blue (dark) +, x, and * symbols indicate results for 
(3 = and e = 0.0125,0.025,0.05 in 2 , respectively. The red 
(bright) circled symbols ©, ®, and © are correspondingly 
for e = —0.0125, —0.025, —0.05 iKi. Each series contains 150 
points. As can be seen, probing with weakly squeezed vac- 
uum can lead to higher purity than probing with unsqueezed 
coherent light. 



environment carries less information on the atomic state. 

Finally, we can also consider very weak squeezing, e = 
±0.0125 iK2- The purity is then above 95% for all d/2 for 
all simulations. For such weak squeezing, however, the 
peaks in the n-distribution are considerably wider than 
for the coherent state scheme. It may be noted that for 
the very weak squeezing, the phase of e apparently has 
a pronounced effect on the purity even for the highest 
values of d/2, while for stronger squeezing, the difference 
disappears as d/2 is increased. 



V. CONCLUSION 

In conclusion, we have analyzed a protocol, which uses 
optical measurements to conditionally prepare an ensem- 
ble of spin- 1/2 atoms in a superposition of a state with 
most of the atoms in the spin up state and a state with 



most of the atoms in the spin down state. The physics 
of the protocol may be understood by referring to the 
plot of the cavity field in Fig. [51 A large horizontal dis- 
tance between two states means that the states are easy 
to distinguish in an x-quadrature measurement, and the 
vertical distance between the states n and — n determines 
the rate of decrease of purity due to light field losses. In- 
creasing 2g/n moves the coherent states around the edge 
of a circle towards the origin, and for an initially coher- 
ent spin state, this increases the purity obtained after 
a given probing time and a favorable measurement out- 
come, but it also decreases the distance d between the 
two peaks obtained in the distribution over the different 
J z eigenstates, because for large 2g/n only states with 
n very close to zero are easily distinguished from states 
with n close to ± J. This is, however, of no concern for 
J = 1, and in this case perfect pure quantum superpo- 
sition states are obtained in the limit 2g/n — > oo. An 
increased number of atoms increases d and the purity of 
the final states, because more states are added to the 
left in Fig. ^ but the increase in d is small, and d/N 
decreases with N. Both the rate of gain of information, 
the rate of loss of photons, and the rate of spontaneous 
emission scale linearly with the probe beam intensity. 

We have also investigated probing with a continu- 
ous beam of squeezed vacuum. This system is signif- 
icantly harder to analyze because the back-action ef- 
fect of the measurements is more complicated, but we 
have found numerically that it is possible to achieve a 
higher purity of the generated states if the probe beam is 
slightly squeezed. The stochastic master equation (jTHJ) is 
quite general and may be used to analyze probing with 
squeezed light in other settings as well. 

The phase shift of the light field per round trip in the 
cavity has been assumed to be infinitesimally small, but 
we note that additional possibilities arise if the phase 
shift per round trip is comparable to 27r. A pure super- 
position of | J) and | — J) can, for instance, be obtained 
if the phase shift per round trip is ±7r for n — =F J and 
light field losses within the atomic ensemble are negligi- 
ble, since, in that case, there is no difference between the 
cavity field strength for the n — ±J atomic states, but 
only a difference between the local phase variation of the 
fields inside the atomic medium. Similarly, if the phase 
shift is 7T per round trip per atom and N is even, it is 
possible to generate pure superpositions of states with 
even values of n and pure superpositions of states with 
odd values of n. 

A coherent spin state is a convenient starting point, 
because it has the desired initial symmetry and is rela- 
tively easy to prepare experimentally, but it could also 
be interesting to consider other initial states. In [25[, it 
has, for instance, been suggested to spin squeeze the co- 
herent spin state in order to increase the weight of the 
larger \n\ components. Another interesting possibility is 
to start from a time-reversal eigenstate. For such a state 
the additional spin symmetries in H are not required, cf. 
Sec. IIVB1 and it is possible to perform arbitrary spin 
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rotations during the probing without breaking the ±n Acknowledgments 
symmetry. This opens the way to use a feedback scheme 
similar to the scheme analyzed in [33| to force the atoms 
into an equal superposition of J and — J (for the case 
where the initial state is an eigenstate of J 2 with eigen- 
value h 2 J(J + 1)). Note that |0)(0| (and any rotated This work was supported by the European Commission 
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